Revealing myopathy spectrum: integrating transcriptional and clinical features of human skeletal muscles with varying health conditions.

Myopathy refers to a large group of heterogeneous, rare muscle diseases. Bulk RNA-sequencing has been utilized for the diagnosis and research of these diseases for many years. However, the existing valuable sequencing data often lack integration and clinical interpretation. In this study, we integrated bulk RNA-sequencing data from 1221 human skeletal muscles (292 with myopathies, 929 controls) from both databases and our local samples. By applying a method similar to single-cell analysis, we revealed a general spectrum of muscle diseases, ranging from healthy to mild disease, moderate muscle wasting, and severe muscle disease. This spectrum was further partly validated in three specific myopathies (97 muscles) through clinical features including trinucleotide repeat expansion, magnetic resonance imaging fat fraction, pathology, and clinical severity scores. This spectrum helped us identify 234 genuinely healthy muscles as unprecedented controls, providing a new perspective for deciphering the hallmark genes and pathways among different myopathies. The newly identified featured genes of general myopathy, inclusion body myositis, and titinopathy were highly expressed in our local muscles, as validated by quantitative polymerase chain reaction.

Haplotype information of large neuromuscular disease genes provided by linked-read sequencing has a potential to increase diagnostic yield.

Rare or novel missense variants in large genes such as TTN and NEB are frequent in the general population, which hampers the interpretation of putative disease-causing biallelic variants in patients with sporadic neuromuscular disorders. Often, when the first initial genetic analysis is performed, the reconstructed haplotype, i.e. phasing information of the variants is missing. Segregation analysis increases the diagnostic turnaround time and is not always possible if samples from family members are lacking. To overcome this difficulty, we investigated how well the linked-read technology succeeded to phase variants in these large genes, and whether it improved the identification of structural variants. Linked-read sequencing data of nemaline myopathy, distal myopathy, and proximal myopathy patients were analyzed for phasing, single nucleotide variants, and structural variants. Variant phasing was successful in the large muscle genes studied. The longest continuous phase blocks were gained using high-quality DNA samples with long DNA fragments. Homozygosity increased the number of phase blocks, especially in exome sequencing samples lacking intronic variation. In our cohort, linked-read sequencing added more information about the structural variation but did not lead to a molecular genetic diagnosis. The linked-read technology can support the clinical diagnosis of neuromuscular and other genetic disorders.

Loss-of-function variants in JPH1 cause congenital myopathy with prominent facial involvement.

Background: Weakness of facial, ocular, and axial muscles is a common clinical presentation in congenital myopathies caused by pathogenic variants in genes encoding triad proteins. Abnormalities in triad structure and function resulting in disturbed excitation-contraction coupling and Ca 2+ homeostasis can contribute to disease pathology. Methods: We analysed exome and genome sequencing data from three unrelated individuals with congenital myopathy characterised by striking facial, ocular, and bulbar involvement. We collected deep phenotypic data from the affected individuals. We analysed the RNA-seq data of one proband and performed gene expression outlier analysis in 129 samples. Results: The three probands had remarkably similar clinical presentation with prominent facial, ocular, and bulbar features. Disease onset was in the neonatal period with hypotonia, poor feeding, cleft palate and talipes. Muscle weakness was generalised but most prominent in the lower limbs with facial weakness also present. All patients had myopathic facies, bilateral ptosis, ophthalmoplegia and fatiguability. While muscle biopsy on light microscopy did not show any obvious morphological abnormalities, ultrastructural analysis showed slightly reduced triads, and structurally abnormal sarcoplasmic reticulum. DNA sequencing identified three unique homozygous loss of function variants in JPH1 , encoding junctophilin-1 in the three families; a stop-gain (c.354C>A; p.Tyr118*) and two frameshift (c.373del p.Asp125Thrfs*30 and c.1738del; p.Leu580Trpfs*16) variants. Muscle RNA-seq showed strong downregulation of JPH1 in the F3 proband. Conclusions: Junctophilin-1 is critical to the formation of skeletal muscle triad junctions by connecting the sarcoplasmic reticulum and T-tubules. Our findings suggest that loss of JPH1 results in a congenital myopathy with prominent facial, bulbar and ocular involvement. Key message: This study identified novel homozygous loss-of-function variants in the JPH1 gene, linking them to a unique form of congenital myopathy characterised by severe facial and ocular symptoms. Our research sheds light on the critical impact on junctophilin-1 function in skeletal muscle triad junction formation and the consequences of its disruption resulting in a myopathic phenotype. What is already known on this topic: Previous studies have shown that pathogenic variants in genes encoding triad proteins lead to various myopathic phenotypes, with clinical presentations often involving muscle weakness and myopathic facies. The triad structure is essential for excitation-contraction (EC) coupling and calcium homeostasis and is a key element in muscle physiology. What this study adds and how this study might affect research practice or policy: This study establishes that homozygous loss-of-function mutations in JPH1 cause a congenital myopathy predominantly affecting facial and ocular muscles. This study also provides clinical insights that may aid the clinicians in diagnosing similar genetically unresolved cases.

Titin copy number variations associated with dominant inherited phenotypes.

BACKGROUND: Titinopathies are caused by mutations in the titin gene (TTN). Titin is the largest known human protein; its gene has the longest coding phase with 364 exons. Titinopathies are very complex neuromuscular pathologies due to the variable age of onset of symptoms, the great diversity of pathological and muscular impairment patterns (cardiac, skeletal muscle or mixed) and both autosomal dominant and recessive modes of transmission. Until now, only few CNVs in TTN have been reported without clear genotype-phenotype associations. METHODS: Our study includes eight families with dominant titinopathies. We performed next-generation sequencing or comparative genomic hybridisation array analyses and found CNVs in the TTN gene. We characterised these CNVs by RNA sequencing (RNAseq) analyses in six patients' muscles and performed genotype-phenotype inheritance association study by combining the clinical and biological data of these eight families. RESULTS: Seven deletion-type CNVs in the TTN gene were identified among these families. Genotype and RNAseq results showed that five deletions do not alter the reading frame and one is out-of-reading frame. The main phenotype identified was distal myopathy associated with contractures. The analysis of morphological, clinical and genetic data and imaging let us draw new genotype-phenotype associations of titinopathies. CONCLUSION: Identifying TTN CNVs will further increase diagnostic sensitivity in these complex neuromuscular pathologies. Our cohort of patients enabled us to identify new deletion-type CNVs in the TTN gene, with unexpected autosomal dominant transmission. This is valuable in establishing new genotype-phenotype associations of titinopathies, mainly distal myopathy in most of the patients.

Novel Compound Heterozygous Splice-Site Variants in TPM3 Revealed by RNA Sequencing in a Patient with an Unusual Form of Nemaline Myopathy: A Case Report.

BACKGROUND: Pathogenic variants in the TPM3 gene, encoding slow skeletal muscle α-tropomyosin account for less than 5% of nemaline myopathy cases. Dominantly inherited or de novo missense variants in TPM3 are more common than recessive loss-of-function variants. The recessive variants reported to date seem to affect either the 5' or the 3' end of the skeletal muscle-specific TPM3 transcript. OBJECTIVES: The aim of the study was to identify the disease-causing gene and variants in a Finnish patient with an unusual form of nemaline myopathy. METHODS: The genetic analyses included Sanger sequencing, whole-exome sequencing, targeted array-CGH, and linked-read whole genome sequencing. RNA sequencing was done on total RNA extracted from cultured myoblasts and myotubes of the patient and controls. TPM3 protein expression was assessed by Western blot analysis. The diagnostic muscle biopsy was analyzed by routine histopathological methods. RESULTS: The patient had poor head control and failure to thrive, but no hypomimia, and his upper limbs were clearly weaker than his lower limbs, features which in combination with the histopathology suggested TPM3-caused nemaline myopathy. Muscle histopathology showed increased fiber size variation and numerous nemaline bodies predominantly in small type 1 fibers. The patient was found to be compound heterozygous for two splice-site variants in intron 1a of TPM3: NM_152263.4:c.117+2_5delTAGG, deleting the donor splice site of intron 1a, and NM_152263.4:c.117 + 164 C>T, which activates an acceptor splice site preceding a non-coding exon in intron 1a. RNA sequencing revealed inclusion of intron 1a and the non-coding exon in the transcripts, resulting in early premature stop codons. Western blot using patient myoblasts revealed markedly reduced levels of the TPM3 protein. CONCLUSIONS: Novel biallelic splice-site variants were shown to markedly reduce TPM3 protein expression. The effects of the variants on splicing were readily revealed by RNA sequencing, demonstrating the power of the method.

Extension of the DNAJB2a isoform in a dominant neuromyopathy family.

Recessive mutations in the DNAJB2 gene, encoding the J-domain co-chaperones DNAJB2a and DNAJB2b, have previously been reported as the genetic cause of progressive peripheral neuropathies, rarely involving pyramidal signs, parkinsonism and myopathy. We describe here a family with the first dominantly acting DNAJB2 mutation resulting in a late-onset neuromyopathy phenotype. The c.832 T > G p.(*278Glyext*83) mutation abolishes the stop codon of the DNAJB2a isoform resulting in a C-terminal extension of the protein, with no direct effect predicted on the DNAJB2b isoform of the protein. Analysis of the muscle biopsy showed reduction of both protein isoforms. In functional studies, the mutant protein mislocalized to the endoplasmic reticulum due to a transmembrane helix in the C-terminal extension. The mutant protein underwent rapid proteasomal degradation and also increased the turnover of co-expressed wild-type DNAJB2a, potentially explaining the reduced protein amount in the patient muscle tissue. In line with this dominant negative effect, both wild-type and mutant DNAJB2a were shown to form polydisperse oligomers.

Mutation update for the ACTN2 gene.

ACTN2 encodes alpha-actinin-2, a protein expressed in human cardiac and skeletal muscle. The protein, located in the sarcomere Z-disk, functions as a link between the anti-parallel actin filaments. This important structural protein also binds N-terminal titins, and thus contributes to sarcomere stability. Previously, ACTN2 mutations have been solely associated with cardiomyopathy, without skeletal muscle disease. Recently, however, ACTN2 mutations have been associated with novel congenital and distal myopathy. Previously reported variants are in varying locations across the gene, but the potential clustering effect of pathogenic locations is not clearly understood. Further, the genotype-phenotype correlations of these variants remain unclear. Here we review the previously reported ACTN2-related molecular and clinical findings and present an additional variant, c.1840-2A>T, that further expands the mutation and phenotypic spectrum. Our results show a growing body of clinical, genetic, and functional evidence, which underlines the central role of ACTN2 in the muscle tissue and myopathy. However, limited segregation and functional data are available to support the pathogenicity of most previously reported missense variants and clear-cut genotype-phenotype correlations are currently only demonstrated for some ACTN2-related myopathies.

Adult-onset dominant muscular dystrophy in Greek families caused by Annexin A11.

OBJECTIVE: Mutations in the prion-like domain of RNA binding proteins cause dysfunctional stress responses and associated aggregate pathology in patients with neurogenic and myopathic phenotypes. Recently, mutations in ANXA11 have been reported in patients with amyotrophic lateral sclerosis and multisystem proteinopathy. Here we studied families with an autosomal dominant muscle disease caused by ANXA11:c.118G > T;p.D40Y. METHODS: We performed deep phenotyping and exome sequencing of patients from four large Greek families, including seven affected individuals with progressive muscle disease but no family history of multi-organ involvement or ALS. RESULTS: In our study, all patients presented with an autosomal dominant muscular dystrophy without any Paget disease of bone nor signs of frontotemporal dementia or Parkinson's disease. Histopathological analysis showed rimmed vacuoles with annexin A11 accumulations. Electron microscopy analysis showed myofibrillar abnormalities with disorganization of the sarcomeric structure and Z-disc dissolution, and subsarcolemmal autophagic material with myeloid formations. Molecular genetic analysis revealed ANXA11:c.118G > T;p.D40Y segregating with the phenotype. INTERPRETATION: Although the pathogenic mechanisms associated with p.D40Y mutation in the prion-like domain of Annexin A11 need to be further clarified, our study provides robust and clear genetic evidence to support the expansion of the phenotypic spectrum of ANXA11.

Comprehensive transcriptomic analysis shows disturbed calcium homeostasis and deregulation of T lymphocyte apoptosis in inclusion body myositis.

OBJECTIVE: Inclusion body myositis (IBM) has an unclear molecular etiology exhibiting both characteristic inflammatory T-cell activity and rimmed-vacuolar degeneration of muscle fibers. Using in-depth gene expression and splicing studies, we aimed at understanding the different components of the molecular pathomechanisms in IBM. METHODS: We performed RNA-seq on RNA extracted from skeletal muscle biopsies of clinically and histopathologically defined IBM (n = 24), tibial muscular dystrophy (n = 6), and histopathologically normal group (n = 9). In a comprehensive transcriptomics analysis, we analyzed the differential gene expression, differential splicing and exon usage, downstream pathway analysis, and the interplay between coding and non-coding RNAs (micro RNAs and long non-coding RNAs). RESULTS: We observe dysregulation of genes involved in calcium homeostasis, particularly affecting the T-cell activity and regulation, causing disturbed Ca2+-induced apoptotic pathways of T cells in IBM muscles. Additionally, LCK/p56, which is an essential gene in regulating the fate of T-cell apoptosis, shows increased expression and altered splicing usage in IBM muscles. INTERPRETATION: Our analysis provides a novel understanding of the molecular mechanisms in IBM by showing a detailed dysregulation of genes involved in calcium homeostasis and its effect on T-cell functioning in IBM muscles. Loss of T-cell regulation is hypothesized to be involved in the consistent observation of no response to immune therapies in IBM patients. Our results show that loss of apoptotic control of cytotoxic T cells could indeed be one component of their abnormal cytolytic activity in IBM muscles.

Bi-allelic loss-of-function OBSCN variants predispose individuals to severe recurrent rhabdomyolysis.

Rhabdomyolysis is the acute breakdown of skeletal myofibres in response to an initiating factor, most commonly toxins and over exertion. A variety of genetic disorders predispose to rhabdomyolysis through different pathogenic mechanisms, particularly in patients with recurrent episodes. However, most cases remain without a genetic diagnosis. Here we present six patients who presented with severe and recurrent rhabdomyolysis, usually with onset in the teenage years; other features included a history of myalgia and muscle cramps. We identified 10 bi-allelic loss-of-function variants in the gene encoding obscurin (OBSCN) predisposing individuals to recurrent rhabdomyolysis. We show reduced expression of OBSCN and loss of obscurin protein in patient muscle. Obscurin is proposed to be involved in sarcoplasmic reticulum function and Ca2+ handling. Patient cultured myoblasts appear more susceptible to starvation as evidenced by a greater decreased in sarcoplasmic reticulum Ca2+ content compared to control myoblasts. This likely reflects a lower efficiency when pumping Ca2+ back into the sarcoplasmic reticulum and/or a decrease in Ca2+ sarcoplasmic reticulum storage ability when metabolism is diminished. OSBCN variants have previously been associated with cardiomyopathies. None of the patients presented with a cardiomyopathy and cardiac examinations were normal in all cases in which cardiac function was assessed. There was also no history of cardiomyopathy in first degree relatives, in particular in any of the carrier parents. This cohort is relatively young, thus follow-up studies and the identification of additional cases with bi-allelic null OBSCN variants will further delineate OBSCN-related disease and the clinical course of disease.

Dominant Distal Myopathy 3 (MPD3) Caused by a Deletion in the HNRNPA1 Gene.

BACKGROUND AND OBJECTIVES: To determine the genetic cause of the disease in the previously reported family with adult-onset autosomal dominant distal myopathy (myopathy, distal, 3; MPD3). METHODS: Continued clinical evaluation including muscle MRI and muscle pathology. A linkage analysis with single nucleotide polymorphism arrays and genome sequencing were used to identify the genetic defect, which was verified by Sanger sequencing. RNA sequencing was used to investigate the transcriptional effects of the identified genetic defect. RESULTS: Small hand muscles (intrinsic, thenar, and hypothenar) were first involved with spread to the lower legs and later proximal muscles. Dystrophic changes with rimmed vacuoles and cytoplasmic inclusions were observed in muscle biopsies at advanced stage. A single nucleotide polymorphism array confirmed the previous microsatellite-based linkage to 8p22-q11 and 12q13-q22. Genome sequencing of three affected family members combined with structural variant calling revealed a small heterozygous deletion of 160 base pairs spanning the second last exon 10 of the heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1) gene, which is in the linked region on chromosome 12. Segregation of the mutation with the disease was confirmed by Sanger sequencing. RNA sequencing showed that the mutant allele produces a shorter mutant mRNA transcript compared with the wild-type allele. Immunofluorescence studies on muscle biopsies revealed small p62 and larger TDP-43 inclusions. DISCUSSION: A small exon 10 deletion in the gene HNRNPA1 was identified as the cause of MPD3 in this family. The new HNRNPA1-related phenotype, upper limb presenting distal myopathy, was thus confirmed, and the family displays the complexities of gene identification.

Out-of-Frame Mutations in ACTN2 Last Exon Cause a Dominant Distal Myopathy With Facial Weakness.

BACKGROUND AND OBJECTIVES: To clinically, genetically, and histopathologically characterize patients presenting with an unusual combination of distal myopathy and facial weakness, without involvement of upper limb or shoulder girdle muscles. METHODS: Two families with a novel form of actininopathy were identified. Patients had been followed up over 10 years. Their molecular genetic diagnosis was not clear after extensive investigations, including analysis of candidate genes and FSHD1-related D4Z4 repeats. RESULTS: Patients shared a similar clinical phenotype and a common pattern of muscle involvement. They presented with a very slowly progressive myopathy involving anterior lower leg and facial muscles. Muscle MRI finding showed complete fat replacement of anterolateral compartment muscles of the lower legs with variable involvement of soleus and gastrocnemius but sparing thigh muscles. Muscle biopsy showed internalized nuclei, myofibrillar disorganization, and rimmed vacuoles. High-throughput sequencing identified in each proband a heterozygous single nucleotide deletion (c.2558del and c.2567del) in the last exon of the ACTN2 gene. The deletions are predicted to lead to a novel but unstructured slightly extended C-terminal amino acid sequence. DISCUSSION: Our findings indicate an unusual form of actininopathy with specific molecular and clinical features. Actininopathy should be considered in the differential diagnosis of distal myopathy combined with facial weakness.

Missense mutations in small muscle protein X-linked (SMPX) cause distal myopathy with protein inclusions.

Using deep phenotyping and high-throughput sequencing, we have identified a novel type of distal myopathy caused by mutations in the Small muscle protein X-linked (SMPX) gene. Four different missense mutations were identified in ten patients from nine families in five different countries, suggesting that this disease could be prevalent in other populations as well. Haplotype analysis of patients with similar ancestry revealed two different founder mutations in Southern Europe and France, indicating that the prevalence in these populations may be higher. In our study all patients presented with highly similar clinical features: adult-onset, usually distal more than proximal limb muscle weakness, slowly progressing over decades with preserved walking. Lower limb muscle imaging showed a characteristic pattern of muscle involvement and fatty degeneration. Histopathological and electron microscopic analysis of patient muscle biopsies revealed myopathic findings with rimmed vacuoles and the presence of sarcoplasmic inclusions, some with amyloid-like characteristics. In silico predictions and subsequent cell culture studies showed that the missense mutations increase aggregation propensity of the SMPX protein. In cell culture studies, overexpressed SMPX localized to stress granules and slowed down their clearance.

Genotype-phenotype correlations in recessive titinopathies.

PURPOSE: High throughput sequencing analysis has facilitated the rapid analysis of the entire titin (TTN) coding sequence. This has resulted in the identification of a growing number of recessive titinopathy patients. The aim of this study was to (1) characterize the causative genetic variants and clinical features of the largest cohort of recessive titinopathy patients reported to date and (2) to evaluate genotype-phenotype correlations in this cohort. METHODS: We analyzed clinical and genetic data in a cohort of patients with biallelic pathogenic or likely pathogenic TTN variants. The cohort included both previously reported cases (100 patients from 81 unrelated families) and unreported cases (23 patients from 20 unrelated families). RESULTS: Overall, 132 causative variants were identified in cohort members. More than half of the cases had hypotonia at birth or muscle weakness and a delayed motor development within the first 12 months of life (congenital myopathy) with causative variants located along the entire gene. The remaining patients had a distal or proximal phenotype and a childhood or later (noncongenital) onset. All noncongenital cases had at least one pathogenic variant in one of the final three TTN exons (362-364). CONCLUSION: Our findings suggest a novel association between the location of nonsense variants and the clinical severity of the disease.

Identification and Characterization of Splicing Defects by Single-Molecule Real-Time Sequencing Technology (PacBio).

Although DNA-sequencing is the most effective procedure to achieve a molecular diagnosis in genetic diseases, complementary RNA analyses are often required.Reverse-Transcription polymerase chain reaction (RT-PCR) is still a valuable option when the clinical phenotype and/or available DNA-test results address the diagnosis toward a gene of interest or when the splicing effect of a single variant needs to be assessed.We use Single-Molecule Real-Time sequencing to detect and characterize splicing defects and single nucleotide variants in well-known disease genes (DMD, NF1, TTN). After proper optimization, the procedure could be used in the diagnostic setting, simplifying the workflow of cDNA analysis.

Is Gene-Size an Issue for the Diagnosis of Skeletal Muscle Disorders?

Human genes have a variable length. Those having a coding sequence of extraordinary length and a high number of exons were almost impossible to sequence using the traditional Sanger-based gene-by-gene approach. High-throughput sequencing has partly overcome the size-related technical issues, enabling a straightforward, rapid and relatively inexpensive analysis of large genes.Several large genes (e.g. TTN, NEB, RYR1, DMD) are recognized as disease-causing in patients with skeletal muscle diseases. However, because of their sheer size, the clinical interpretation of variants in these genes is probably the most challenging aspect of the high-throughput genetic investigation in the field of skeletal muscle diseases.The main aim of this review is to discuss the technical and interpretative issues related to the diagnostic investigation of large genes and to reflect upon the current state of the art and the future advancements in the field.

Improved Criteria for the Classification of Titin Variants in Inherited Skeletal Myopathies.

BACKGROUND: Extensive genetic screening results in the identification of thousands of rare variants that are difficult to interpret. Because of its sheer size, rare variants in the titin gene (TTN) are detected frequently in any individual. Unambiguous interpretation of molecular findings is almost impossible in many patients with myopathies or cardiomyopathies. OBJECTIVE: To refine the current classification framework for TTN-associated skeletal muscle disorders and standardize the interpretation of TTN variants. METHODS: We used the guidelines issued by the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) to re-analyze TTN genetic findings from our patient cohort. RESULTS: We identified in the classification guidelines three rules that are not applicable to titin-related skeletal muscle disorders; six rules that require disease-/gene-specific adjustments and four rules requiring quantitative thresholds for a proper use. In three cases, the rule strength need to be modified. CONCLUSIONS: We suggest adjustments are made to the guidelines. We provide frequency thresholds to facilitate filtering of candidate causative variants and guidance for the use and interpretation of functional data and co-segregation evidence. We expect that the variant classification framework for TTN-related skeletal muscle disorders will be further improved along with a better understanding of these diseases.

Recurrent TTN metatranscript-only c.39974-11T>G splice variant associated with autosomal recessive arthrogryposis multiplex congenita and myopathy.

We present eight families with arthrogryposis multiplex congenita and myopathy bearing a TTN intron 213 extended splice-site variant (NM_001267550.1:c.39974-11T>G), inherited in trans with a second pathogenic TTN variant. Muscle-derived RNA studies of three individuals confirmed mis-splicing induced by the c.39974-11T>G variant; in-frame exon 214 skipping or use of a cryptic 3' splice-site effecting a frameshift. Confounding interpretation of pathogenicity is the absence of exons 213-217 within the described skeletal muscle TTN N2A isoform. However, RNA-sequencing from 365 adult human gastrocnemius samples revealed that 56% specimens predominantly include exons 213-217 in TTN transcripts (inclusion rate ≥66%). Further, RNA-sequencing of five fetal muscle samples confirmed that 4/5 specimens predominantly include exons 213-217 (fifth sample inclusion rate 57%). Contractures improved significantly with age for four individuals, which may be linked to decreased expression of pathogenic fetal transcripts. Our study extends emerging evidence supporting a vital developmental role for TTN isoforms containing metatranscript-only exons.